% BEM-LAP-MAT project
% Matlab/Freemat codes
% Copyright 2008 Stephen Kirkup 
% http://www.researchgate.net/profile/Stephen_Kirkup
% University of Central Lancashire
% Issued under the GNU General Public License 2007, see gpl.txt
% www.boundary-element-method.com
% contact: stephen@boundary-element-method.com
% http://www.researchgate.net/profile/Stephen_Kirkup


% function [vertpts,elemvert]=square 
% Returns a list of edges of the panels that make up the square
%  in vertpts and the list of elements that make up the boundary
%  are defined by the indices of the two vertices in elemvert.
% Note that the vertices of each element in elemvert must be
%  defined in the clockwise direction.

function [vertpts,elemvert]=square

  vertpts=[0.00000000000 0.00000000000;
           0.00000000000 0.01250000000;
           0.00000000000 0.02500000000;
           0.00000000000 0.03750000000;
           0.00000000000 0.05000000000;
           0.00000000000 0.06250000000;
           0.00000000000 0.07500000000;
           0.00000000000 0.08750000000;

           0.00000000000 0.10000000000;
           0.01250000000 0.10000000000;
           0.02500000000 0.10000000000;
           0.03750000000 0.10000000000;
           0.05000000000 0.10000000000;
           0.06250000000 0.10000000000;
           0.07500000000 0.10000000000;
           0.08750000000 0.10000000000;

           0.10000000000 0.10000000000;
           0.10000000000 0.08750000000;
           0.10000000000 0.07500000000;
           0.10000000000 0.06250000000;
           0.10000000000 0.05000000000;
           0.10000000000 0.03750000000;
           0.10000000000 0.02500000000;
           0.10000000000 0.01250000000;

           0.10000000000 0.00000000000; 
           0.08750000000 0.00000000000;
           0.07500000000 0.00000000000;
           0.06250000000 0.00000000000;
           0.05000000000 0.00000000000;
           0.03750000000 0.00000000000;
           0.02500000000 0.00000000000;
           0.01250000000 0.00000000000];
           
           elemvert=[1:32; 2:33]
           elemvert(2,32)=1